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Abstract 

A  model  for  a  nickel-metal  hydride  cell  was  constructed  based  on  the  planar  electrode  approximation.  The  mass  balances  of  active 
species,  the  kinetics  of  electrochemical  reactions,  the  ohmic  effects  of  internal  resistance,  and  the  energy  balance  of  the  whole  cell  were 
considered  in  the  model.  An  empirical  approach  was  utilized  to  account  for  the  hysteresis  potential  behavior  of  the  nickel  electrode.  The 
model  predictions  showed  favorable  agreement  with  the  experimental  data.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

The  nickel-metal  hydride  battery  is  one  of  the  latest 
battery  technologies.  The  nickel-metal  hydride  battery 
has  some  advantages,  e.g.  higher  energy  density  and  envir¬ 
onmental  friendliness,  over  the  more  commonly  used 
rechargeable  batteries  (such  as  the  lead-acid  battery  and 
the  nickel-cadmium  battery).  However,  it  also  has  some 
disadvantages,  e.g.  higher  cost  and  higher  self-discharge. 

A  nickel-metal  hydride  cell  module  (Fig.  1)  consists  of 
three  major  components;  the  positive  nickel  electrode,  the 
nylon  separator,  and  the  negative  metal  hydride  electrode. 
Generally,  both  the  nickel  electrode  and  the  metal  hydride 
electrode  are  thin  porous  electrodes.  In  a  nickel-metal 
hydride  cell,  the  nickel  electrode  determines  the  cell  capa¬ 
city  and  the  metal  hydride  electrode  determines  the  cell 
cycle  life. 

Normally,  the  metal  hydride  electrode  is  pre-charged  (i.e. 
has  residue  charge  when  the  cell  is  fully  discharged)  to  avoid 
oxygen  gas  generation  during  over  discharge  and  has  extra 
capacity  to  avoid  hydrogen  gas  generation  during  over¬ 
charge.  Since  the  metal  hydride  material  gradually  loses 
capacity  through  usage  due  to  the  corrosion  by  the  KOH 
electrolyte,  the  extra  capacity  in  the  metal  hydride  electrode 
also  reduces  the  rate  at  which  the  capacity  is  lost  in  the  cell. 

The  electrolyte  in  a  nickel-metal  hydride  cell  is  concen¬ 
trated  KOH  solution,  which  has  good  electric  conductivity 
for  a  wide  range  of  temperatures.  Some  LiOH  is  usually 
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added  to  the  electrolyte  to  improve  the  performance  of  the 
electrodes. 

In  the  normal  operation  of  a  nickel-metal  hydride  cell, 
there  are  four  major  reactions,  as  listed  below. 

Positive  nickel  electrode; 

discharge 

NiOOH  +  H20  +  e~  ?±  Ni(OH)2  +  OH  (1) 

charge 

20H~  -4  1 02  +  H20  +  2e~  (2) 

Negative  metal  hydride  electrode: 

discharge 

MH  +  OH'  H20  +  M  +  e~  (3) 

charge 

\  02  +  H20  +  2e~  — *  20H~  (4) 

There  is  a  main  and  a  side  reaction  at  each  of  the  electrodes. 
The  main  reaction  at  the  positive  electrode  is  the  redox 
reaction  of  the  nickel  active  material  and  at  the  negative 
electrode  the  redox  reaction  of  the  metal  hydride  material. 
The  side  reaction  at  the  positive  electrode  is  oxygen  evolu¬ 
tion  and  at  the  negative  electrode  oxygen  reduction.  The  two 
side  reactions  are  coupled  together  by  the  oxygen  transport 
from  the  positive  electrode  to  the  negative  electrode. 

With  the  rapid  development  of  computer  technology, 
modeling  and  simulation  has  become  an  indispensable 
approach  in  the  design  and  optimization  of  many  engineer¬ 
ing  systems.  For  batteries,  the  modeling  and  simulation  can 
be  used  to:  (i)  analyze  the  influence  of  some  physical 
parameters  (e.g.  electrode  thickness  and  particle  size) 
on  battery  behavior  for  design  optimizations  and  scale-up 
investigations;  (ii)  predict  battery  behavior  under  different 
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Nomenclature 

dnQg 

specific  surface  area  of  the  negative  electrode 
(cm2/cm3) 

^pos 

specific  surface  area  of  the  positive  electrode 
(cm2/cm3) 

^cell 

total  surface  area  of  the  cell  (cm2) 

A 

^neg 

total  geometry  area  of  the  negative  electrode 
(cm2) 

A 

^pos 

total  geometry  area  of  the  positive  electrode 
(cm2) 

concentration  of  KOH  electrolyte  (mol/cm3) 

^e,ref 

reference  concentration  of  KOH  electrolyte 
(mol/cm3) 

CH+ 

concentration  of  nickel  hydroxide  in  nickel 
active  material  (mol/cm3) 

^H+  ,max 

maximum  concentration  of  nickel  hydroxide 
in  nickel  active  material  (mol/cm3) 

cH+,ref 

reference  concentration  of  nickel  hydroxide  in 
nickel  active  material  (mol/cm3) 

CMH 

concentration  of  hydrogen  in  metal  hydride 
material  (mol/cm3) 

^MH,max 

maximum  concentration  of  hydrogen  in  metal 
hydride  material  (mol/cm3) 

cMH,ref 

reference  concentration  of  hydrogen  in  metal 
hydride  material  (mol/cm3) 

CP 

heat  capacity  of  the  cell  (J/g/K) 

Ea,k 

activation  energy  of  reaction  k  (J) 

f 

function  of  reactant  concentrations 

F 

Faraday’s  constant  (96,487  C/eq) 

h 

heat  transfer  coefficient  between  the  cell  and 
the  environment  (W/cm2/K) 

^cell 

current  of  the  nickel-metal  hydride  cell  (A) 

io  ,k 

exchange  current  density  of  reaction  k  (A/cm2) 

lo,k,ref 

exchange  current  density  of  reaction  k  at 
reference  reactant  concentrations  (A/cm2) 

^k,ref 

limiting-current  density  of  reaction  k  at 
reference  reactant  concentrations  (A/cm2) 

jk 

current  density  of  reaction  k  (A/cm2) 

kc 

parameter  in  the  empirical  hysteresis  equation 
for  the  charge  process 

kd 

parameter  in  the  empirical  hysteresis  equation 
for  the  discharge  process 

kQ 

parameter  in  the  empirical  hysteresis  equation 
for  the  open-circuit  process 

^neg 

thickness  of  the  negative  electrode  (cm) 

^pos 

thickness  of  the  positive  electrode  (cm) 

ly,  neg 

effective  thickness  of  the  layer  of  metal 
hydride  material  (cm) 

^y,pos 

effective  thickness  of  the  layer  of  nickel  active 
material  (cm) 

^MH 

loading  of  metal  hydride  material  (g/cm2) 

-^Ni(OH)2 

loading  of  nickel  active  material  (g/cm2) 

^cell 

mass  of  the  cell  (g) 

Po2 

oxygen  pressure  in  a  nickel-metal  hydride  cell 
(atm) 

PO2  ,ref 

reference  oxygen  pressure  in  a  nickel-metal 
hydride  cell  (atm) 

R 

universal  gas  constant  (8.3143  J/mol/K  or 
82.0562  atm  cm3/mol/K) 

^int 

internal  ohmic  resistance  of  the  nickel-metal 
hydride  cell  (Q) 

soc 

state  of  charge 

t 

time  (s) 

T 

temperature  (K) 

To 

room  temperature  (298.15  K) 

uk 

equilibrium  potential  of  reaction  k  at  standard 
conditions  (V) 

f/l,c 

apparent  open-circuit  potential  of  the  redox 
reaction  of  nickel  active  material  at  standard 
conditions  during  the  whole  range  charge 
process  (V) 

f/l,d 

apparent  open-circuit  potential  of  the  redox 
reaction  of  nickel  active  material  at  standard 
conditions  during  the  whole  range  discharge 
process  (V) 

^cell 

potential  of  a  nickel-metal  hydride  cell  (V) 

y 

y  gas 

gas  volume  in  a  nickel-metal  hydride  cell  (cm3) 

Greeks  letters 

^0neg 

potential  difference  at  the  solid-liquid  inter¬ 
face  on  the  negative  electrode  (V) 

pos 

potential  difference  at  the  solid-liquid  inter¬ 
face  on  the  positive  electrode  (V) 

tfq 

apparent  open-circuit  potential  at  the  starting 
point  of  a  process  (V) 

0eq  ,k 

equilibrium  potential  of  reaction  k  (V) 

0eq,l,c 

apparent  open-circuit  potential  of  the  redox 
reaction  of  nickel  active  material  during  the 
whole  range  charge  process  (V) 

0eq,l,d 

apparent  open-circuit  potential  of  the  redox 
reaction  of  nickel  active  material  during  the 
whole  range  discharge  process  (V) 

4*h,c 

transient  hysteresis  potential  of  the  redox 
reaction  of  nickel  active  material  during  the 
charge  process  (V) 

4*h,d 

transient  hysteresis  potential  of  the  redox 
reaction  of  nickel  active  material  during  the 
discharge  process  (V) 

4*  h,o 

transient  hysteresis  potential  of  the  redox 
reaction  of  nickel  active  material  during  the 
open-circuit  process  (V) 

0neg,l 

potential  in  the  liquid  electrolyte  on  the 
negative  electrode  (V) 

0neg,s 

potential  in  the  solid  metal  hydride  material  on 
the  negative  electrode  (V) 

^pos,l 

potential  in  the  liquid  electrolyte  on  the 
positive  electrode  (V) 

0  pos,s 

potential  in  the  solid  nickel  active  material  on 
the  positive  electrode  (V) 

PNi(OH)2 

density  of  nickel  active  material  (g/cm3) 

P  MH 

density  of  metal  hydride  material  (g/cm3) 
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Fig.  1.  Schematic  diagram  of  a  NiMH  cell  module. 


operating  conditions  (e.g.  charge/discharge  rate  and  cut-off 
voltage)  for  operation  optimizations;  and  (iii)  investigate 
unclear  physical/chemical  mechanisms  by  comparing  the 
model  predictions  with  extensive  experimental  data. 

The  nickel  electrode  has  been  used  and  researched  for  a 
long  time  [1-6].  However,  well-established  theories  for 
some  characteristics  of  the  nickel  electrode  (e.g.  second 
discharge  plateau  [7]  and  hysteresis  potential  behavior  [8]) 
are  still  not  currently  available.  A  significant  amount  of 
modeling  work  has  been  done  over  the  years  on  the  nickel 
electrode  either  alone  or  in  a  battery  [9-14],  The  most 
comprehensive  model  for  the  nickel  electrode  is  a 
pseudo-2D  model  with  considerations  of  both  the  processes 
across  the  macroscopic  electrode  dimension  and  the  pro¬ 
cesses  inside  the  microscopic  nickel  active  material  layer. 

Some  modeling  efforts  on  the  metal  hydride  electrode 
have  been  conducted  in  recent  years  [15-18].  A  pseudo-2D 
model  of  the  metal  hydride  electrode  similar  to  that  of  the 
nickel  electrode  was  also  utilized. 

However,  very  limited  modeling  work  is  reported  in  the 
literature  for  the  complete  nickel-metal  hydride  cell.  Paxton 
and  Newman  [19]  developed  a  discharge  model  of  a  nickel- 
metal  hydride  cell.  They  found  that  the  proton  diffusion  in 
the  nickel  active  material  could  be  neglected  at  normal 
charge/discharge  rates.  Based  on  their  model,  some  optimi¬ 
zations  of  cell  design  with  respect  to  the  electrode  thickness 
and  porosity  were  conducted.  Later,  Gu  et  al.  [20]  con¬ 
structed  a  more  comprehensive  model  for  a  nickel-metal 
hydride  cell.  The  oxygen  reactions  on  both  the  positive  and 
negative  electrodes  were  included  in  their  model.  Using  their 
model,  the  effects  of  the  oxygen  evolution  rate  and  the 
charge/discharge  rate  on  the  cell  behavior  were  investigated. 

The  model  developed  by  Paxton  and  Newman  [19]  and 
the  one  by  Gu  et  al.  [20]  are  both  isothermal  ones.  However, 
in  reality,  significant  heat  effects  are  involved  in  the  opera¬ 
tion  of  a  battery,  and  many  processes  in  a  battery  are  very 
sensitive  to  temperature  changes.  Thus,  it  is  highly  desirable 
to  include  the  thermal  effects  in  a  nickel-metal  hydride  cell 
model. 

Furthermore,  these  two  models  used  a  reversible  electro¬ 
chemical  reaction  to  describe  the  behavior  of  the  nickel 
active  material.  However,  the  potential  behavior  of  the 
nickel  active  material  shows  significant  hysteresis  features 


that  cannot  be  described  by  a  reversible  electrochemical 
reaction.  Thus,  for  some  operating  conditions,  there  may  be 
significant  prediction  errors  with  these  two  models. 

In  general,  the  number  of  details  to  be  considered  in  a 
nickel-metal  hydride  cell  will  determine  the  complexity  of 
the  model  to  be  constructed.  The  more  details  considered  in 
the  model,  the  more  parameters  and  equations  will  appear, 
and  the  more  complex  and  computational  intensive  the 
model  program  will  become.  For  the  nickel-metal  hydride 
battery,  the  following  hierarchy  of  models  could  be  con¬ 
structed: 

1 .  Empirical  model:  fitted  equations  from  experimental  data. 

2.  Lumped  model  (with  the  planar  electrode  approxima¬ 
tion):  equations  of  mass  balance,  charge  balance, 
electrochemical  kinetics,  and  energy  balance,  etc. 

3.  Distributed  model  (with  ID  or  pseudo-2D  porous 
electrode  models):  equations  in  (2)  plus  equations  of 
distributed  processes. 

4.  Coupled  electrochemical/thermal  model:  equations  in 
(1),  (2),  and  (3)  plus  the  heat  transfer  equations  in 
multiple  dimensions. 

The  primary  objective  of  this  work  is  to  provide  a  model 
of  a  nickel-metal  hydride  cell  for  power  system  simulations. 
Since,  many  models  of  different  devices  need  to  be  solved 
together  in  a  system  simulation,  the  battery  model  must  be 
simple  enough  to  make  the  whole  system  simulation  com¬ 
putationally  feasible.  As  a  result,  the  lumped  modeling 
approach  was  used  in  this  work.  The  advantages  of  this 
approach  are  (i)  the  model  is  based  on  first-principle 
mechanisms;  (ii)  model  equations  are  simple  to  solve; 
and  (iii)  the  model  can  be  easily  incorporated  into  a  system 
model  or  a  thermal  model.  The  disadvantages  of  this 
approach  are  (i)  model  assumptions  may  not  be  valid  under 
some  situations;  and  (ii)  effects  of  many  cell  parameters 
cannot  be  evaluated. 


2.  Governing  equations  of  the  model 


The  model  presented  here  uses  the  planar  electrode 
approximation,  i.e.  each  electrode  is  treated  as  a  huge  planar 
electrode  with  a  surface  area  equal  to  the  total  internal 
surface  area  of  the  porous  electrode.  Since,  battery  electro¬ 
des  are  normally  made  very  thin  to  minimize  ohmic  effects,  a 
planar  electrode  model  usually  gives  insignificant  errors  and 
costs  much  less  computation  compared  to  a  rigorous  porous 
electrode  model. 

The  kinetics  of  reactions  (Eqs.  ( 1)— (3))  are  described  with 
the  Butler- Volmer  equation. 


—  exp 


0.5  F 
RT 


(A^pos 


0eq,l) 


(5) 
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72  =  io,2 


73  —  l*o,3 


eXp((^(A</Vg  -  Yq.3^ 

-  exP  (A0neg  -  <^eq,3)^) 


(6) 


(7) 


The  potential  differences  at  the  solid-liquid  interface  on  the 
electrodes  are  given  by 

Af/lp„s  l/lpos,s  ^/ipos.l  (^1 


A</>neg  =  ^neg.s  -  </>neg,l 


(9) 


The  exchange  current  densities  in  Eqs.  (5)-(7)  can  be 
calculated  from  the  reference  states 


1 0,1  —  !o,l, ref 


{  CH+ 

VCH+,, 


ref 


C’e.rcf 


°.5  /  X  0.5 

'  I  [ '  .max  '  H 


/£a.l  (\  1 

exphrlr-7b 

2 

/  C-e  ' 

^o,2  —  ^o,2, ref  1 

V^e.ref 


PO2 

02,ref 


^H+,max  ^H+,ref. 


’  (E&  ( 1  1 

expU(r-7b 


^o,3  —  ^o,3,  ref 


(  CMH 


7 


^MH,max  ^MH 


\£-MH,ref/  \^e,ref/  \OvlH,max  ^MH,ref 

/£a.3  /  1  1 

expUlr-7o 


(10) 


(11) 

0.5 


(12) 


For  the  nickel  active  material,  the  hysteresis  potential  beha¬ 
vior  exists,  i.e.  the  charge  and  discharge  potential  curves 
form  hysteresis  loops.  In  this  work,  an  empirical  approach 
has  been  used  to  characterize  the  hysteresis  potential.  For  the 
whole  range  of  charge  and  discharge  processes,  two  appar¬ 
ent  open-circuit  potential  curves  based  on  the  Nernst 
equation  are  utilized. 


_JT  ,(T  ,d£uc  RT 
q.i.c  —  UU  +  (T  —  7o)  I  —  In 


cH+ 


cH+ 


ceCH+ 


(13) 


_  jj  .  frp  rj-,  X  dt/14  ,  RT  ) 

q.t.d  —  17i,d  +  (T  —  T0)  —  +  —  In 


CH+ 


CH+ 


ceCH+ 


(14) 


However,  when  switching  between  different  processes,  there 
is  a  relaxation  process  before  the  ideal  charge/discharge 
potential  curves  are  reached,  as  shown  in  Figs.  2  and  3.  Three 
empirical  equations  are  used  in  this  work  to  represent  this 
relaxation  process. 

K  =  Cq  +  (<7q,l,c  ^  Cq)(l  "  (15) 


Fig.  2.  Qualitative  diagram  of  hysteresis  potential  behavior  for  the  charge 
process. 


Ka  =  <Yq  +  (</>eq,l,d  -  Cq)(l  "  e  1  ^ 1  ( ^  ^  >  ) 

(16) 

^h,o  =  Cq  +  (^eq,l,d-^q)(l-e-^)) 

(17) 

In  the  above  equations  <p°q  is  the  value  at  the  starting  time. 
The  values  for  parameters  kc ,  kd,  and  kQ  can  be  determined 

from  the  experimental  data. 

The  equilibrium  potential  for  the  oxygen 
given  by 

reactions  is 

,  „  ,(T  T  x  dU2  RT ,  (Po?\ 

K,2  -  t72  +  (r  -  r0)  —  +  —  In  I 

(18) 

while  the  equilibrium  potential  for  the  metal  hydride  reac¬ 
tion  is  expressed  by  Paxton  and  Newman  [19] 


^eq,3  =  u3  +  (T  -  To)  ^  ^  lnce  +  9.712  x  IQ-4 


0.23724  exp  - 


28.057cmh\ 


Cmh  ,max 

2.7302  x  10"4 


J 


[(cmh/ cmh  ,max)  H"~  0.010768] 


(19) 


soc 

Fig.  3.  Qualitative  diagram  of  hysteresis  potential  behavior  for  the 
discharge  process. 


B.  Wu  et  al.  /  Journal  of  Power  Sources  101  (2001)  149-157 


153 


Due  to  the  huge  electrochemical  driving  force  for  oxygen 
reduction  on  the  negative  electrode,  a  limiting-current 
equation  is  used  for  the  rate  of  reaction  Eq.  (4). 


P02 

J  4  = - 

PO2  ,ref 


i-l.rei  CXp 


(20) 


The  charge  balances  on  the  electrodes  are  described  by 


fed  I  — A  pos^pos^pos  (/  \  T  J2)  (21) 

fcell  —  ^neg^neg^neg  C/3  T  7*4  )  (22) 

The  mass  balance  of  the  nickel  active  material  is  given  by 


dci 


H+ 


7t 
'  F 


(23) 


^,pos  dr 

where  the  effective  thickness  of  the  nickel  active  material  is 


d4i(OH)2 


Pn{(  OH)  2  (pos  ^pos 


(24) 


Substitution  of  Eq.  (33)  into  Eq.  (30)  yields 
dT  , 

CpWlceU  —  ^^celll,-*  -^a)  ^'cell^cell 

<34> 

*=1,4  x  / 

The  SOC  of  the  cell  is  given  by  the  charged  state  of  the 
nickel  active  material 

SOC  =  1  -  Ch+  (35) 

CH+  ,max 

Given  the  applied  current  on  the  cell  /cen,  there  are  seven 
dependent  variables  (vceii,  cH+,  cmh,  po2,  Ar/>pos,  A(/>neg,  and 
T)  that  need  to  be  solved  with  the  model.  The  procedure  used 
for  serial  simulation  of  electrochemical  systems,  based  on 
two  FORTRAN  solvers  DAEIS  and  DASRT,  has  been  used 
to  simulate  the  nickel-metal  hydride  cell  undergoing  many 
processes  [21]. 


The  mass  balance  of  metal  hydride  material  is  given  by 


Ocmi 
;  dr 


74 

'  F 


(25) 


where  the  effective  thickness  of  metal  hydride  active 
material  is 


(26) 


(27) 


,  _  Lmh 

(v.neg  —  ~ 

PMH‘negflneg 

The  mass  balance  of  oxygen  is  described  by 

Eg  as  dpo2  (4  pos^posfpog/2  I  AnegGEneg/neg/4) 

~RT~dt~ ~  F 

The  separator  is  treated  as  an  ohmic  resistance,  which  gives 

r'cell^int  =  0pos,l  —  ^neg.l  (28) 

The  cell  potential  is  calculated  by 

Cl-]]  0 pos.s  —  ^neg.s  A(j) po^  ^fineg  3”  Cell^int  (29) 

The  energy  balance  of  the  whole  cell  is  described  by 
dT  , 

<7>mceii  —  hAcen(T  Ta)  4-  /'ccib’ccii 


-  J2>k 


Jeq,k 


-  T 


d  0, 


eq,& 


dr 


Since,  the  following  equations  hold 

Kk=Uk  +  (T~To)^  +  Tf(0 

d<^eq,*  d  Uk  ,  v 

~sT=  dr+/(c) 

we  have 


(30) 

(31) 

(32) 


3.  Results  and  discussion 

The  parameter  values  used  in  the  model  are  those  of  a 
4.5  Ah  (rated  capacity)  nickel-metal  hydride  cell  produced 
by  VARTA  Battery  AG  (www.varta.com).  These  values  are 
listed  in  Table  1 .  The  actual  cell  capacity  after  some  cycling 
is  close  to  5.3  Ah.  The  experimental  data  was  generated 
during  testing  at  the  Ford  Motor  Company  Scientific 
Research  Laboratory,  in  Dearborn,  Michigan,  USA,  using 
a  BTS-600  Firing  Circuits  battery  tester  and  software. 
Where  possible,  the  model  predicted  values  and  the  actual 
test  values  are  overlain  in  the  figures  presented. 

The  cell  potential  with  respect  to  time  for  three  processes 
is  shown  in  Fig.  4.  The  high  potential  plateau  at  the  end  of 
the  charge  process  corresponds  to  the  overcharge  period 
when  the  oxygen  evolution  and  reduction  reactions  dom¬ 
inate.  Since,  the  rate  of  oxygen  evolution  is  much  slower 


d(/), 


W  T  dT 


eq  ,k 


uk 


To 


dUk 
d  T 


(33) 


Fig.  4.  Cell  potential  for  the  following  processes  (model  predictions  in 
line):  0.8  A  charge  for  8  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to 
1.0  V. 
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Table  1 

Parameters  used  for  the  model  of  a  nickel-metal  hydride  cell 


Parameter 

Value 

flpos  (cm2/cm3) 

4000.0 

flneg  (cm2/cm3) 

3000.0 

A>os  (cm  ) 

325.0 

ncg  (cm-) 

360.0 

Aceii  (cm2) 

77.0 

ce  (mol/cm3) 

7.0  X  l(r3 

Ceref  (mol/cnT) 

1.0  x  10~3 

fMH.raax  (mol/cm3) 

1.0  x  10  1 

CMH.ref  (mol/cm3) 

0-5  (CMH,max) 

cH+  max  (mol/cm3) 

3.7  x  10~2 

CH+,ref 

0-5  (CNi(OH)2,max) 

to.Uef  (A/cm2) 

1.0  x  10~4 

(o,2.ref  (A/cm-) 

2.0  x  10-" 

<0,3,ref  (A/cm2) 

1.0  x  10~4 

*4, ref  (A/cm2) 

1.0  x  10~4 

kc 

1.0  x  10~3 

kd 

1.0  x  10~3 

ko 

1.0  x  10~2 

Ea  i  (kJ/mol) 

10 

Z?a.2  (kJ/mol) 

120 

£a  3  (kJ/mol) 

10 

Ea.4  (kJ/mol) 

10 

Ipos  (cm) 

3.3  x  10~2 

(lie g  (cm) 

2.8  x  10~2 

^Ni(OH),  (g/cm2) 

6.8  x  10~2 

Lmh  (g/cm2) 

1.13  x  10~‘ 

PO2 ,ref  (atm) 

1.0 

Kin,  (fi) 

5.0  x  10~3 

To  (K) 

298.15 

(A, c  (V) 

0.527 

t/,,d  (V) 

0.427 

u2  (V) 

0.4011 

£/3(V) 

-0.8279 

£4(V) 

0.4011 

dUi/dT  (mV/K) 

-1.35 

dU2/dT  (mV/K) 

-1.68 

dt/3/dr  (mV/K) 

-1.55 

dUJdT  (mV/K) 

-1.68 

Vgas  (cm3) 

1.0  x  10~' 

PNi(OH),  (g/cm3) 

3.4 

Pmh  (g/cm3) 

7.49 

than  that  of  the  nickel  reaction,  a  higher  overpotential  is 
needed  to  maintain  the  same  charge  rate.  It  is  observed  that 
the  cell  was  overcharged  for  more  than  2  h.  The  full  capacity 
of  the  cell  can  be  easily  determined  from  the  discharge 
process  of  this  operation. 

The  cell  temperature  with  respect  to  time  is  shown  in 
Fig.  5.  The  heat  generation  at  the  early  stage  of  the  charge 
process  is  insignificant.  During  the  overcharge  period,  since 
almost  all  the  input  energy  is  wasted  on  heat  generation, 
there  is  a  rapid  increase  in  cell  temperature.  At  the  end  of  the 
discharge  process,  the  depletion  of  reactants  results  in  high 
overpotential  loss  that  also  causes  a  cell  temperature 
increase.  It  is  worth  mentioning  that  hydrogen  absorption 
into  the  metal  hydride  material  is  an  exothermal  process, 
while  hydrogen  release  from  the  metal  hydride  material  is  an 


Fig.  5.  Cell  temperature  for  the  following  processes  (model  predictions  in 
line):  0.8  A  charge  for  8  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to 
1.0  V. 


endothermic  process.  This  has  significant  influence  on  the 
thermal  behavior  of  the  whole  nickel-metal  hydride  cell  and 
can  be  observed  at  low  charge/discharge  rates. 

The  mole  fractions  of  the  charged  species  on  the  positive 
and  negative  electrodes  with  respect  to  time  are  shown  in 
Fig.  6.  During  the  early  period  of  the  charge  process  and  the 
entire  discharge  process,  nearly  constant-slope  lines  are 
obtained  for  these  mole  fractions.  This  indicates  that  side 
reactions  are  insignificant  during  these  periods. 

A  5%  pre-charge  on  the  metal  hydride  electrode  is  used  in 
the  simulation.  It  should  be  noted  that  less  than  half  of  the 
metal  hydride  material  is  utilized.  During  the  overcharge 
period,  there  is  little  if  any  increase  in  cell  capacity.  This  is 
due  to  the  dominance  of  oxygen  side  reactions  that  consume 
almost  all  the  input  current. 

The  cell  potential  with  respect  to  time  for  three  different 
operations  is  shown  in  Fig.  7.  In  these  three  operations,  the 
open-circuit  and  discharge  processes  are  the  same,  but  the 
charge  rate  is  different  for  all  three.  While  there  is  good 


Fig.  6.  Model  predicted  mole  fractions  of  charged  active  material  on  the 
positive  and  negative  electrodes  for  the  following  processes:  0.8  A  charge 
for  8  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to  1.0  V. 
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Time  (h) 

Fig.  7.  Cell  potential  for  three  operations  (model  predictions  in  line):  (i) 
20.0  A  charge  for  0.22  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to 
1.0  V;  (ii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  4.0  A 
discharge  to  1.0  V;  and  (iii)  0.8  A  charge  for  5.5  h,  open-circuit  for  0.5  h, 
and  4.0  A  discharge  to  1.0  V. 


correlation  between  the  model  predictions  and  the  experi¬ 
mental  data  for  the  charge  and  discharge  processes,  some 
errors  are  noted  for  the  open-circuit  process.  This  indicates 
that  efforts  to  improve  the  model  should  concentrate  on 
providing  a  better  description  of  the  open-circuit  process. 

At  high  charge  rates,  due  to  the  high  ohmic  and  over- 
potential  losses,  high  charge  potential  is  needed  as  expected. 
It  is  worth  mentioning  that  the  Varta  nickel-metal  hydride 
cell  performs  well  even  at  a  30  C  charge  rate.  Thus,  the  solid 
phase  diffusion  limitation  is  unlikely  to  be  important  in  the 
cell  at  normal  charge  rates. 

The  cell  temperature  with  respect  to  time  for  the  same 
operations  is  given  in  Fig.  8.  At  high  charge  rates,  the  cell 
heat  generation  is  significant  due  to  high  ohmic  and  over- 
potential  losses,  which  may  cause  the  cell  temperature  rise 
more  than  10°C.  Since,  the  cell  is  not  fully  charged  for  the 


Fig.  8.  Cell  temperature  for  three  operations  (model  predictions  in  line): 
(i)  20.0  A  charge  for  0.22  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to 
1.0  V;  (ii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  4.0  A 
discharge  to  1.0  V;  and  (iii)  0.8  A  charge  for  5.5  h,  open-circuit  for  0.5  h, 
and  4.0  A  discharge  to  1 .0  V. 


Time  (h) 

Fig.  9.  Model  predicted  oxygen  pressure  in  the  cell  for  three  operations: 
(i)  20.0  A  charge  for  0.22  h,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to 
1.0  V;  (ii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  4.0  A 
discharge  to  1.0  V;  and  (iii)  0.8  A  charge  for  5.5  h,  open-circuit  for  0.5  h, 
and  4.0  A  discharge  to  1 .0  V. 


given  operations,  it  can  be  expected  that  overcharge  will 
cause  a  even  more  significant  rise  of  the  cell  temperature  at 
high  charge  rates,  which  may  cause  thermal  safety  problems. 

The  model  predicted  oxygen  pressure  in  the  cell  is  given 
in  Fig.  9.  It  is  observed  that  at  high  charge  rates,  the  oxygen 
pressure  rise  can  be  significant  at  the  end  of  the  charge 
process.  There  is  negligible  oxygen  pressure  during  the 
discharge  process. 

The  cell  potential  with  respect  to  time  for  three  different 
operations  is  shown  in  Fig.  10.  In  these  operations  the 
charge,  and  open-circuit  processes  are  the  same,  but  the 
discharge  process  is  different  for  all  three  operations.  As 
expected,  low  discharge  potential  is  resulted  at  high  dis¬ 
charge  rates  due  to  the  high  ohmic  and  overpotential  losses. 
A  S-bend  discharge  curve  is  shown  at  all  three  discharge 
rates.  Since,  the  cell  is  not  discharged  from  the  fully  charged 


Fig.  10.  Cell  potential  for  three  operations  (model  predictions  in  line):  (i) 
4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  20.0  A  discharge  to 
1 .0  V;  (ii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  4.0  A 
discharge  to  1 .0  V;  and  (iii)  4.0  A  charge  for  1 . 1  h,  open-circuit  for  0.5  h, 
and  0.8  A  discharge  to  1.0  V. 
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Fig.  11.  Cell  temperature  for  three  operations  (model  predictions  in  line): 
(i)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  20.0  A  discharge  to 
1.0  V;  (ii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h,  and  4.0  A 
discharge  to  1.0  V;  and  (iii)  4.0  A  charge  for  1.1  h,  open-circuit  for  0.5  h, 
and  0.8  A  discharge  to  1.0  V. 


state,  the  Nernst  equation  will  predict  a  lower  potential  at  the 
early  discharge  period.  The  higher  potential  at  the  early 
stage  of  the  discharge  process  is  actually  an  effect  of  the 
hysteresis  phenomenon,  i.e.  potential  relaxation  when 
switching  processes. 

The  cell  temperature  for  the  same  operations  is  given  in 
Fig.  11.  At  high  discharge  rates,  the  cell  heat  generation  is 
significant  due  to  high  ohmic  and  overpotential  losses,  and 
that  causes  sharp  cell  temperature  rises.  Thus,  good  thermal 
management  is  needed  if  high  discharge  rates  are  used  to 
avoid  a  thermal  run-away  situation. 

The  cell  potential  with  respect  to  time  for  a  series  of  pulse 
charge/open-circuit  operations  followed  by  an  open-circuit 
and  a  discharge  process  is  shown  in  Fig.  12.  The  cell 
potential  with  respect  to  time  for  a  charge  process  followed 
by  an  open-circuit,  and  then  a  series  of  pulse  discharge/open- 
circuit  operations  is  shown  in  Fig.  13.  The  potential  differ¬ 
ence  between  the  pulse  charge/discharge  processes  and  the 


Fig.  12.  Cell  potential  for  the  following  processes  (model  predictions  in 
thick  line):  repeat  (4.0  A  charge  for  0.1 1  h  and  open-circuit  for  0.167  h)  for 
11  times,  open-circuit  for  0.5  h,  and  4.0  A  discharge  to  1.0  V. 


Fig.  13.  Cell  potential  for  the  following  processes  (model  predictions  in 
thick  line):  4.0  A  charge  for  1.1  h  and  open-circuit  for  0.5  h,  and  repeat 
(4.0  A  discharge  for  0.11  h  and  open-circuit  for  0.167  h)  for  10  times,  and 
4.0  A  discharge  to  1 .0  V. 


following  open-circuit  processes  is  mainly  the  effects  of 
ohmic  resistance.  However,  careful  comparison  shows  that 
the  open-circuit  potential  is  different  between  pulse  charge 
and  pulse  discharge  operations  at  corresponding  cell  SOC. 
This  potential  difference  gives  a  clear  demonstration  of  the 
significance  of  the  hysteresis  phenomenon.  This  behavior 
has  been  satisfactorily  reflected  with  the  model. 

The  cell  potential  with  respect  to  time  for  five  different 
operations  is  shown  in  Fig.  14.  In  these  operations,  the  cell  is 
charged  to  different  SOC  levels  (at  the  same  charge  rate)  and 
then  it  is  discharged  at  the  same  rate  after  the  same  open- 
circuit  period.  It  is  observed  that  the  cell  discharge  potential 
is  higher  than  that  predicted  by  the  Nernst  equation,  espe¬ 
cially  when  discharging  from  a  lower  SOC.  This  is  a  typical 


Fig.  14.  Cell  potential  for  the  following  five  different  operations  (model 
predictions  in  thick  line):  (i)  4.0  A  charge  for  14  min,  open-circuit  for 
30  min,  and  4.0  A  discharge  to  1.0  V;  (ii)  4.0  A  charge  for  22  min,  open- 
circuit  for  30  min,  and  4.0  A  discharge  to  1.0  V;  (iii)  4.0  A  charge  for 
38  min,  open-circuit  for  30  min,  and  4.0  A  discharge  to  1.0  V;  (iv)  4.0  A 
charge  for  46  min,  open-circuit  for  30  min,  and  4.0  A  discharge  to  1.0  V; 
and  (v)  4.0  A  charge  for  54  min,  open-circuit  for  30  min,  and  4.0  A 
discharge  to  1 .0  V. 
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Fig.  15.  Cell  potential  for  the  following  processes  (model  predictions  in 
thick  line):  4.0  A  charge  for  80  min,  open-circuit  for  30  min,  4.0  A 
discharge  for  8  min,  open-circuit  for  30  min,  4.0  A  charge  for  8  min,  open- 
circuit  for  30  min,  4.0  A  discharge  for  16  min,  open-circuit  for  30  min, 
4.0  A  charge  for  16  min,  open-circuit  for  30  min,  4.0  A  discharge  for 
24  min,  open-circuit  for  30  min,  4.0  A  charge  for  24  min,  open-circuit  for 
30  min,  4.0  A  discharge  for  32  min,  open-circuit  for  30  min,  4.0  A  charge 
for  32  min,  open-circuit  for  30  min,  4.0  A  discharge  for  40  min,  open- 
circuit  for  30  min,  4.0  A  charge  for  40  min,  open-circuit  for  30  min,  4.0  A 
discharge  for  48  min,  open-circuit  for  30  min,  4.0  A  charge  for  48  min, 
open-circuit  for  30  min,  4.0  A  discharge  for  56  min,  open-circuit  for 
30  min,  4.0  A  charge  for  56  min,  open-circuit  for  30  min,  4.0  A  discharge 
for  64  min,  open-circuit  for  30  min,  4.0  A  charge  for  64  min,  open-circuit 
for  30  min,  4.0  A  discharge  to  1.0  V  cut-off  voltage. 


effect  of  hysteresis  phenomenon.  Though  some  prediction 
errors  exist,  the  model  predictions  give  reasonable  qualita¬ 
tive  trends. 

The  cell  potential  with  respect  to  time  for  many  processes 
is  shown  in  Fig.  15.  The  cell  is  fully  charged  first.  Then  the 
following  operations  are  repeated  on  the  cell:  open-circuit 
for  half  an  hour,  discharge  partially,  open-circuit  for  half  an 
hour,  and  fully  charge  again.  The  depth  of  discharge  gra¬ 
dually  increases  for  each  subsequent  operation.  It  is  found 
that  the  potential  of  the  charge  process  is  lower  than  that 
predicted  by  the  Nernst  equation,  another  typical  effect  of 
the  hysteresis  phenomenon.  However,  there  is  good  correla¬ 
tion  between  the  model  predictions  and  the  experimental 
data  for  these  operations. 

The  processes  presented  in  Figs.  12-15  are  very  common 
operations  for  a  nickel-metal  hydride  cell.  However,  pre¬ 
dicting  the  corresponding  potential  behavior  of  the  cell  is 
very  difficult  due  to  the  influence  of  the  hysteresis  phenom¬ 
enon.  Consequently,  a  model  that  uses  a  reversible  nickel 
reaction,  will  encounter  significant  prediction  errors.  The 
model  developed  in  this  work,  using  the  empirical  approach, 
reflects  satisfactorily  the  hysteresis  potential  behavior  of  the 
nickel-metal  hydride  cell. 


4.  Conclusions 

A  lumped  model  of  a  nickel-metal  hydride  cell  was 
developed  with  attention  given  to  many  important  mechan¬ 
isms  in  the  cell.  One  important  characteristic  of  the  nickel 
electrode,  the  hysteresis  potential  behavior,  is  described  with 
empirical  equations.  The  model  can  be  used  to  predict  the 
potential,  pressure,  and  temperature  behavior  of  a  nickel- 
metal  hydride  cell  for  charge/open-circuit/discharge  pro¬ 
cesses.  The  model  predictions  correlated  well  with  the 
experimental  data  of  a  VARTA  nickel-metal  hydride  cell. 
Further  research  work  is  desirable  to  provide  a  theoretical 
explanation  of  the  hysteresis  phenomenon  of  the  nickel 
electrode. 
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